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Abstract True biaxial tests of granular materials are 
investigated by applying the principle of minimal dissi- 
pation and comparing to two dimensional contact dy- 
namics simulations. It is shown that the macroscopic 
steady state manifested by constant stress ratio and 
constant volume is the result of the ever changing mi- 
croscopic structure which minimizes the dissipation rate. 
The shear band angle in the varying shear band struc- 
tures is found to be constant. We also show that in- 
troducing friction on the walls reduces the degeneracy 
of the optimal shear band structures to one for a wide 
range of parameters which gives a non-constant stress 
ratio curve with varying aspect ratio that can be calcu- 
lated. 

Keywords biaxial shear, granular flow, shear band, 
least dissipation, optimization, Mohr-Coulomb theory 



1 Introduction 

Shear bands seem inevitable whenever granular mate- 
rial is subject to quasistatic deformation. The form and 
the position of these failure zones depend on the exper- 
imental setup [T] and sometimes even on the specific 
sample [2]. 

Classical continuum theories predict strain fields from 
the stresses of the material under load. Generally ac- 
cepted constitutive relations are still missing for dense 
granular materials, not to mention a microscopic un- 
derstanding of the parameters. 

Therefore it is legitimate to approach the physics of 
shear bands from a radically simplified perspective in 
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which one idealizes shear bands as strain field singulari- 
ties and regards them as the basic dynamical quantities 
of the theory. This is the approach used in this paper 
where shear bands are dividing the granular material 
into domains which are rigidly displaced with respect 
to each other. The realized shear band configuration is 
determined using the principle of least dissipation [31 

HJ. 

We apply this method to the true biaxial tester [SJ 
,6j containing a non-cohesive granular material, where 
steady state flow can be observed. It allows the defor- 
mation of a granular medium in a brick shaped con- 
tainer with perpendicular, frictionless walls. This keeps 
the principal axes of the stress and strain tensors con- 
stant and aligned perpendicular to the container walls. 
We would like to find out what is the structure of the 
fiow that allows for an apparent macroscopic steady 
state, and under what conditions can one predict the 
velocity field of the sample. 

In this paper we present first some numerical simu- 
lation results on true biaxial tests to motivate the the- 
oretical model, which is then solved for the setup with 
frictionless and later with frictional walls. Finally we 
test the predictions on the simulation results. 



2 Simulation results 

2.1 Definition of the system 

In this paper we use distinct element simulations to 
study the shearing of frictional granular material in the 
so called true biaxial tester , which consists of two 
fixed walls normal to the z direction, stress controlled 
confining walls normal to the x direction (normal stress 
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Fig. 1 Snapshots from simulated biaxial tests. Time is progressing from (a) to (c). The particles are represented with arrows 
proportional to their actual velocity. The color (online) is given according to which quadrant the vector is pointing. Courtesy 
of D. Kadau [5]. 



axx)^ and strain controlled compressing walls perpen- 
dicular to the y direction (normal strain rate eyy<0). 

The simulations were done using the contact dy- 
namics method 7] . Since the behavior is invariant along 
the z direction we simulate only a two-dimensional sys- 
tem in the xy plane. Samples of 10000 frictional discs 
were enclosed in a rectangle by four walls. For the con- 
stant confining stress axx on the vertical and compress- 
ing strain rate iyy on the horizontal walls we chose val- 
ues such that the inertial number [S] was always of the 
order of 10"'* [S] so that the deformation can be re- 
garded as quasistatic. Further details of the simulations 
can be found in |S]. The particles were frictional with 
microscopic friction coefficient /iniic=0.3. Tests with dif- 
ferent particle numbers and friction coefficients show 
qualitatively identical results, and thus are not pre- 
sented in this paper. 

First, we present results with frictionless walls, cor- 
responding to the usual experimental setup [B]. 



2.2 Shear band patterns 

For frictionless walls, after an initial transient, the stress 
ayy measured on the strain controlled and the strain 
rate Sxx on the stress controlled walls fluctuate around 
a constant value This was interpreted as a steady 
state, however, the underlying microscopic behavior re- 
mains unknown. 

On Fig. [T] we show some instantaneous velocity snap- 
shots of a simulation at different instances. Visual ob- 



servation of the velocity snapshots clearly indicate the 
existence of shear bands which are narrow straight lines 
reflected by the boundary or absorbed by the corners. 
In a different, two dimensional setup Hall et al. also 
reported complicated shear band structures |10j . 

Simulations with different initial configurations show 
in general different shear band structures at given in- 
stances, however certain easily recognizable configura- 
tions appear around a given aspect ratio and never too 
far from it. 

The above observations give rise to two important 
questions which we want to answer in this paper: (i) 
Why does the stress ratio (Jyy/<7xx remain constant in 
spite of the strongly varying velocity field, (ii) Is it pos- 
sible to predict the shear band patterns in biaxial tests? 
In the next section we try to answer these questions 
based on the principle of minimal dissipation. 



3 Principle of minimal dissipation 

3.1 Theory 

Before introducing the principle of minimal dissipation 
we summarize our assumptions: The system is described 
in the frame, where the center of mass is fixed at the 
origin. We assume that the sample has already reached 
steady state flow |9] and is quasistatic, i.e. there are no 
accelerations. Ignoring the small volume fluctuations al- 
lows the calculation of the velocities of all walls from 
the given strain rate iyy: The upper and lower walls 
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move with velocity — Uj,, respectively Uy, the left and 
right ones with —u^, respectively Ux, where 

Uy ^ Ly\eyy\/2 ^11(1 U ^ ^ L ^ £ X X / '2 ^ L ^l^yy \ / 2 . (l) 

We model the shear bands as infinitely narrow, piece- 
wise straight lines running all the way from one con- 
tainer wall to another without splitting or merging. 
These fines cut the system into domains inside which all 
grain velocities are the same. The domain boundaries 
are either shear band segments or walls. 

The velocity of the particles in a domain must ful- 
fill the following boundary conditions: (i) At a wall the 
normal component of a domain velocity must coincide 
with the velocity of the wall. The parallel component of 
the domain velocity can be freely chosen as the walls are 
frictionless. (ii) A shear band segment requires a discon- 
tinuity of the tangential velocity with some magnitude 
Av, while the normal component must be continuous. 

An immediate consequence of these two conditions 
is that shear bands cannot end at a wall, but must be 
"reflected" under some angle or end in one of the cor- 
ners of the rectangular container. Thus, in the following 
we consider only shear band configurations represented 
by one or more piecewise straight paths B that close 
upon themselves or begin and end in a corner. Despite 
these restrictions the set of admissible shear band con- 
figurations is very large. They will be analyzed in detail 
in the next sections. 

Among these possible configurations we select the 
ones which minimize the energy dissipation rate 



E= ^lAvPdl, 



where the integral extends over the whole path B repre- 
senting the shear band configuration, Av is the tangen- 
tial velocity difference across the shear band, P is the 
normal stress acting on the line element dl. We suppose 
that the material is homogeneous and is described by 
a single effective friction coefficient /i which is constant 
in the system. 



3.2 Diagonal shear band 

We start our analysis with the the simplest possible case 
which is a single shear band spanning diagonally the 
sample [Figs. [TJc) andUJa)]. Let us denote the aspect 
ratio by Ly/Lx=taiif5, the angle between the major 
principal axis (the y-axis in the present case) and the 
normal of the shear band by a (cf. Fig. [2]). In this special 
case a=l3. 

The velocity in the two domains should match the 
bounding walls, so the velocity difference between the 
two domains is Av=Uy/ sina. The unit length on the 




Fig. 2 (a) The simplest shear band configuration, the diag- 
onal, (b) Shear band configuration for a=/3 with one degree 
of freedom. 



shear band is dl—dx/ cos a. P{a) can be derived from 
force balance on the shear band as there are no ac- 
celerations: the horizontal components of the normal 
and tangential force components on the line element are 
P{a)dl sin a and —fj,P{a)dl cos a, respectively. Their sum 
must be equal to Uxxdy, where dy /dl= sina, yielding 



P{a) 



(3) 



1 — /X cot a 

An immediate consequence is that the shear band 
angle a must be larger than the internal friction angle 
(f) = arctan/i, as the normal stress cannot be tensile. 
When a approaches 4> the normal stress P diverges, 
expressing the impossibility that such a system moves 
(at finite confining stress axx)- 

Thus the dissipation rate of the diagonal shear band 
takes the following form: 



(2) Mo^) = 



fJ^liyyWxxV 



COS a sin a(l — fi cot a) 



tan a -I- 1 



= H\Syy\axxV "^" - ' - J-^-j 

tana - /i 

where V—LxLy denotes the volume of the sample. Min- 
imizing the above equation yields 

amin = "MC 4 + '2 ' ' 

where aMC is the Mohr-Coulomb angle. 

The equation ^ is an important result. It says that 
the dissipation rate for the diagonal shear band is min- 
imal when the angle of the shear band equals to the 
Mohr-Coulomb angle. 

The dissipation rate can be identified with the work 
per unit time done by the walls: 



E — ^ {^xx^xx ~t- ^yy^yy) 

From this the stress ratio can be derived: 
E 



(6) 



s{a) = 



1 



1 



cos a sin a(l — /x cot a) 



(7) 
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Fig. 3 Combined shear band configuration samples (a) 
qjr = 2/1 (b) qjr = 3/2. The arrows show the velocities 
of the particles in the domains. 



The stress ratio of this diagonal configuration varies 
with the aspect ratio. Thus the diagonal shear band 
should be observed only around /3 ~ q;mc otherwise 
the stress ratio cannot be constant. 



3.3 Complex shear band structures 

Let us first examine the configuration shown on Fig.l^Tb). 
For this configuration the formula for the dissipation 
rate is the same as in Eq. since the length of the 
shear band is twice but the velocity jump is half than 
for the diagonal shear band. Thus again the minimum 
of the dissipation rate is at a=/3=aMC • 

Let us note, that Fig. [H^b) has a degree of freedom. 
The reflection point on the top wall can be anywhere. 
Moving the reflection point R into the corner reduces 
this configuration to the diagonal shear band. Fig.[2fa) 
is a special case of Fig. ^h). 

Let us call shear band configuration Fig. [^fb) an 
elementary cell. More complicated shear band configu- 
rations can be constructed if we place q times r elemen- 
tary cells in a rectangular lattice by matching the refiec- 
tion point of the shear bands and scaling the velocity 
components to match at the cell boundaries (examples 
are shown in Fig. [S]). The minimization presented in 
[TT] shows that for these configurations the dissipation 
rate is minimal if all shear band angles are equal and 
given by a—auc- Furthermore the value of the dissi- 
pation rate is the same as for a single elementary cell 
(velocity jump and length scale inversely proportion- 
ally), namely: 



(8) 



It is important to stress that the minimal dissipation 
rate for the complex shear band structures is indepen- 
dent of the values of q and r. 

On the other hand if q/rj^l then the aspect ratio of 
the sample is different from tan aMC- Thus for all aspect 
ratios tan /? and for arbitrary small e we can create a 
shear band configuration with jg/r — tanaMc/ tan/3|<e 



which have all the same -Emin with shear band angle 
equaling the Mohr-Coulomb angle. 

The stress ratio is a linear function of the minimal 
dissipation rate cf. Eq. ([7]). Thus in the scope of this 
idealized theory at all aspect ratios we should find a 
different shear band configuration, and at the same time 
observe a constant stress ratio. 

This constant stress ratio can be expressed with the 
effective friction coefficient /i cf. Eqs. (IS|) and ([7|): 

fs-l\ 

LL = tan arcsin (9) 

\s + lj 

The above formula was already derived in a different 
way in [I] . 

The stress ratio can be measured in the simulations 
to be 5=1.54. From this value the effective friction co- 
efficient gets /i=0.213. 

This result answers question (i) of Sec. 12.21 The 
apparent macroscopic steady state is the result of the 
complete lack of microscopic steady state which at each 
aspect ratio favors a different set of shear band config- 
urations with the same minimal dissipation rate and 
stress ratio. 

The answer for question (ii) is also clear within 
the scope of these idealized shear band structures: For 
all rational tan aMC / tan /? there are infinitely many 
pairs of q and r that fulfill tanaMc/tan/J^g/r, and 
are equivalent in terms of dissipation rate. Thus it is 
impossible to predict which shear band configuration 
will be realized. 

We should also note here [11] that the linear na- 
ture of the system allows for linear combinations of 
shear band configurations which also enlarges the set 
of available configurations for a given aspect ratio. An 
example: Fig. [1] (b) is a linear combination of Fig. [T] (c) 
and its mirror image. 

The results shown here are valid for homogeneous 
materials. In reality the effective friction coefficient fluc- 
tuates around its mean and does so the dissipation rate. 
The actual minimum does not only depend on the av- 
erage value but also on the variance. The integral in 
Eq. (HI averages the fluctuations of the effective friction 
coefficient along the path. The overall variance of the 
path is inversly proportional to the length of the shear 
band. This way the shorter shear bands have larger vari- 
ance and thus more likely to be minimal. A consequence 
is that one can easily find nice, simple shear band con- 
figurations like the examples on Fig. [T] in almost all 
simulations. 

In summary in a perfectly homogeneous material 
one would find different shear band configurations at 
each time step. The apparent macroscopic steady state 
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is the result of the fact that for the optimal configu- 
ration the dissipation rate depends only on the shear 
band angle and is independent of the aspect ratio. 



4 Numerical results with frictionless walls 

4.1 Shear band detection 

An important prediction of our theory is that the shear 
band angle is independent of the sample aspect ratio 
and equals to the Mohr-Coulomb angle auc- In order 
to test this result on the simulations we need to detect 
the position of the shear band. Many quantities were 
already introduced to measure the local shear intensity 
[12]. We found that for our purpose the best choice is 
the Frobenius norm of the gradient of the coarse grained 
velocity v (see [l3]), which is defined as follows: 



1000 



iVvl 



id^V,,f + {d^Vy)^ + {dyVcc)^ + idyVy)HW) 



The strain rate tensor can be decomposed into vol- 
umetric £^ and shear e'^ strain part: 



2 



VV , 



-yy 



(11) 



which leads to a decomposition of 



(12) 



where A± = ±\J e^J' + e%y are the eigenvalues of the 
shear strain tensor. 

In this paper we call I^g the local shear intensity. 



4.2 Shear band direction 

In general the angular correlation of I^g is a good way 
to measure the direction of the shear bands marked 
with high values of I^g- However, in our case both a 
and n—a angle shear bands are possible that cross each 
other many times. This adds a lot of noise to the angular 
correlations. Thus we chose another method: 

In order to measure the directions of the shear bands 
a segment 1 of length l=\l\ larger than the width of the 
shear band is placed at point x. Then the correlation 
between the local shear intensity is calculated at the 
end points of the segment only if all along the segment 
1 the valus of I^g are larger than its average lyg. Thus 
the new angular correlation function is defined as: 
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Fig. 4 The segment angular correlation function of the 
shear intensity averaged for 4 samples and 400 timesteps, for 
2.6>tan/3>0.5. The length of the segment was |1|=10 (x), 
resp. |1| = 12 (+) in units of average particle diameter. The 
continuous curve is a Gaussian fit to the latter data for points 
C/>500. 



This definition has the advantage that it avoids contri- 
butions from points lying in different shear band seg- 
ments. 

On Fig. |3]we show the values of Cj for two different 
segment length. It displays a strong peak at an angle of 
aMC=51.5°±l°. The predicted value from the measured 
stress ratio is Q!mc=51°- 

As a side result, the distribution functions predict a 
shear band width of about 6.5 particle diameters, which 
is consistent with the picture of narrow shear bands. 



5 Frictional walls 

5.1 The X shaped shear band structure 

The negative answer to the question (ii) of Sec. 
raises another problem of importance. Namely: Is there 
a case where the principle of minimal dissipation can 
predict the observed shear band configuration? The an- 
swer can only be affirmative, if we can reduce the de- 
generacy of the possible solutions. 

The simplest way to do this is to introduce friction 
on the walls. Wall friction introduces a finite difference 
in the dissipation rate for shear band configurations 
with or without velocity slip at the walls. Furthermore, 
it is obvious that there is only one shear band config- 
uration with no tangential velocity at the walls, since 
the following two restrictions must be met: (i) it may 
not have any shear band reflected at the walls, which 
would mean different tangential velocities in the two do- 
mains, (ii) domains touching different walls should be 
separated by shear bands. These two conditions imply 
that perpendicular velocities in the domains touching 
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Fig. 5 In the top figure the data points show the stress ratio 
as function of the angle of the diagonal. The data points are 
ensemble averages over 30 different realizations. The fitted 
curve is Eq. ([7]) with /x = 0.226. The snapshots presented be- 
low were taken at the indicated instances. The blue color is 
proportional to the shear intensity I^g. Particles are color 
coded according to the quarter they are pointing to: red- 
north, yellow-east, green-south, black-west. 



the walls are only possible for the X-shaped shear band 
configuration, which is composed of the two diagonal 
shear bands, see Fig.[Uc). 

We expect thus that the stress ratio is not constant 
any more but has a minimum at aspect ratio of the 
Mohr-Coulomb angle (tan (3= tan aMc) which is the op- 
timal for the X-shaped shear band structure: Further- 
more close to this point the stress ratio should be de- 
scribed by Eq. ©. 

The above equation can be tested on our data. Since 
iyy is a simulation parameter, the only fit parameter 
we are left with is the effective friction coefficient of the 
bulk fi. The value n = 0.226, Eq. © fits very well the 
numerical data over a wide range of aspect ratios as 
shown on Fig. [5] Moreover snapshots of the system also 
clearly indicate the X-shaped shear band structure. 

The difference between the friction coefficients in 
case of frictionless and frictional walls will be the sub- 
ject of a future paper. 

The answer to the question raised at the beginning 
of this section is thus yes, indeed there are cases when 
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Fig. 6 The forces acting on the domains in case of the di- 
agonal shear band. 



the principle of minimal dissipation is capable of pre- 
dicting the actually observed shear bands. In the next 
section we discuss if we can predict the overall stress 
ratio curve and the other observed shear band configu- 
rations. 



5.2 Complex shear band structures 

The calculation of the dissipation rate for a general con- 
figuration in case of frictional walls is more complicated 
because pressure acting on shear band segments does 
not only depend on the shear band angle. Instead it 
must be calculated from the force balance of the do- 
mains. On Fig. iniwe show an example using the diago- 
nal shear band. 

The force balance shown on Fig. |6] has one external 
parameter, the normal force Nx—UxxLy acting on the 
confining walls, two material parameters: the effective 
friction coefficient of the bulk ^ and the wall friction 
/i^, and two unknowns, the normal force on the shear 
band N^b and the normal force on the compressing walls 
Ny=ayyLx. Thc force balance of the two domains gives 
four (two components for each domain) linear equations 
for the unknowns, but for symmetry reasons only two 
are independent. The solution of these linear equations 
allows the calculation of the dissipation rate and the 
stress ratio. 

The above process sounds simple, but in general it 
has many problems. First, in all cases there are more 
equations than unknowns and it may happen that a 
given shear band configuration often seen in the fric- 
tionless simulations gets realizable only at a given as- 
pect ratio. Second, the analytic solution of the linear 
system of equations is of course possible, but the general 
analytic minimization is impossible due to the trigono- 
metric functions. So we are left with analyzing numer- 
ically the candidate configurations one by one. 
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Fig. 7 Top: A typical shear band structure beyond the X 
shape limit. Bottom: The schematic representation of the in- 
complete shear band structure. 



Instead of testing this infinite set of configurations 
- motivated by the typical shear band structure shown 
on Fig. [7] - we define here an incomplete shear band 
configuration in wliich only part of the shear band sys- 
tem is explicitely defined but still the minimization can 
be carried out. Here we focus on the aspect ratio range 
of /3 ^ ayic where the X shaped shear band becomes 
unfavourable. 

Our assumptions are the following: (i) In order to 
minimize slip on the walls the generic configuration 
must have perpendicular velocities on the largest pos- 
sible surface: This gives an X shape in the middle and 
two half X on the side walls. Since we are in the aspect 
ratio regime tan /3<1, the half X shear bands on the con- 
fining walls can cover the whole wall with an optimal 
angle ai=aMC- The middle X does not cover the whole 
wall and it has an angle of a2 which does not necessar- 
ily match the Mohr- Coulomb angle, (ii) The structure 
in the place left out (gray parts on Fig. [7]) remains un- 
known. We know however, that the force balance is kept 
on the known shear band segments, (iii) The velocity 
slip on the horizontal walls in the gray section is also 
unknown. Here we assume a simple linear dependence 
which gives an average slip velocity oiux/2. Finally, we 
note here, that the visual observations fully support the 
above assumptions as shown on Fig. [T) 

The minimization with respect to a2 of the dissipa- 
tion rate for each aspect ratio is numerically possible if 
we know the effective wall friction coefficient. The lat- 
ter parameter can be determined e.g. from simple shear 
simulations [T3] and it was found to be /iu,~0.12. The 
result of nmnerical minimization at different aspect ra- 
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Fig. 8 Test of the incomplete X configuration on the nu- 
merical data. 



tios match prefectly the simulation data as shown on 
Fig. El 

Thus in conclusion the introduction of friction on 
the walls destroys the apparent macroscopic steady state 
but fixes the shear band structure for a large range of 
aspect ratios. The prediction of the stress ratio is still 
possible with the help of two parameters, the effective 
friction of the bulk and the wall. Beside the central 
regime we cannot predict the exact shear band struc- 
ture just some features of it. 



6 Conclusion 

In conclusion, in this paper we applied the principle 
of minimal dissipation to the biaxial shear test. From 
a theory based on simple assumptions, in case of fric- 
tionless walls, we were able to show that the apparent 
macroscopic constant stress ratio is a surprising result 
of the permanent minimization of the dissipation rate. 
This is manifested by a changing shear band configura- 
tion with constant shear band angle which is equal to 
the Mohr-Coulomb angle. For each sample aspect ratios 
different shear band configurations are optimal, but the 
high degeneracy of these optimal configurations denies 
us the possibility to predict the observed shear band 
configurations. 

The above degeneracy can be reduced to a single op- 
timal configuration by introducing friction on the walls, 
in which case for a wide range of aspect ratios a single 
configuration is optimal which can be clearly observed 
in simulations. Knowing the optimal configuration, or 
at least the important characteristics of it allows us to 
predict the the macroscopic stress ratio. 
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